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Abstract 

We investigate radioactivity from the decay sequence of ^^Ti in young supernova 
remnants (SNRs), Cassiopeia A (Cas A) and SN 1987A. It is shown by a hnear 
analysis that ionization of "^^Ti, a pure electron capture decay isotope, affects the 
radioactivity contradistinctively in these two SNRs: Ionization of ^"^Ti to H-like and 
He-like states enhances its present radioactivity in Cas A, while such high-ionization 
decreases its radioactivity in SN 1987A. We briefly discuss the enhancement factor of 
the present radioactivity of Cas A considering microscopic (atomic/nuclear) physics 
combined with a hydrodynamical SNR evolution model. For SN 1987A, we have 
obtained the initial '^^Ti mass of (0.82 - 2.3) X 1O~^M0 from our Monte-Carlo 
simulations. The resulting fluxes of 7 and hard X-rays emerged from the ^^Ti decay 
are given for the current and future experiments. 
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1 Introduction 



The initial yield of "^^Ti that is synthesized by a single event of a core-collapse 
supernova explosion is very crucial to constrain dynamics of core-collapse su- 
pernova nucleosyntehsis. This is because ^^Ti is synthesized at the vicinity of 
so-called mass cut, that divides the matter which accretes on a compact object 
and ejecta which is scattered into interstellar space. For this, the initial mass 
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of ^Ti depends sensitively on 1) the location of the mass cut, 2) the maxi- 
mum temperature and the maximum density behind the shock wave, and 3) 
the internal structure (< 2 Mq from the center) of a progenitor. 

For the above reason, it is very interesting to compare theoretical predictions 
of the ''^Ti yield with "observed" values. Since ^^Ti is radioactive, we can 
detect its radioactivity and derive the initial '^^Ti mass from it; ^^Ti decays by 
electron capture to "^"^Sc, emitting 67.9 keV and 78.4 keV nuclear deexcitation 
hnes. Then ^Sc decays almost exclusively by positron emission to ^Ca, which 
emits 1.16 MeV deexcitation line. The emitted positron ends up with 511 keV 
annihilation line. So far, detection of the 1.16 MeV line from Cas A with 
COMPTEL/CGRO experiment (e.g., lyudin et al. 1994; Schonfelder et al. 
2000) and reconfirmation of this by 67.9 and 78.4 keV lines with BeppoSAX 
(Vink et al. 2001) have allowed us to do such a comparison. It is also expected 
to detect the ^^Ti nuclear hnes from other young galactic SNRs and SN 1987A 
in LMC in near future. 

The halflife of ^"^Ti is a key quantity for its radioactivity to be an important 
observable in young SNRs, and hence has been intensively studied in labora- 
tories after the first detection of the nuclear 7-ray fiux by lyudin et al. (1994). 
Compilation of recent 8 experiments which were performed after 1998 (see., 
e.g., Hashimoto et al. 2001 and references therein; Fiilop et al. 2000) gives 
weighted mean halflife of ti/2 = 60± 1 yr (the error is 1 a, statistical). It is 
this timescale of the halflife that makes ^^Ti a useful diagnostic isotope. 

However, a crucial point here is that ^^Ti decays only by orbital electron 
capture. This is because the decay Q-value from the ground state of "^^Ti to 
the second excited state of ^"^Sc (branching ratio of 99.3%) is less than twice 
the electron rest mass, which is at least required for positron emission to 
be allowed by producing two 511 keV 7-photons when a positron annihilates 
with an electron (and so does that to the first excited state of "^^Sc for the 
rest of the minor fraction of the branch). Thus we should be careful to apply 
the experimental halflife to this problem, because halflife measurements in 
laboratories are done for neutral atoms: The electric environment for *^^Ti in 
a young SNR may be very much different from that in laboratories. 

As we shall point out later, there is a clear possibility of ionization of ^"^Ti 
ongoing in SN 1987A. Also, there are indications that '^'^Ti in Cas A is highly 
ionized; this may be conflrmed directly in the future spectroscopic observations 
in X-rays. In this article, therefore, we are going to discuss the radioactivity 
of ^^Ti in young SNRs taking the role of ionization into consideration. Previ- 
ous studies are found in Mochizuki et al. (1999) and Mochizuki (2001). In the 
following, a linear analysis is presented in section 2, and a result of the radioac- 
tivity calculated for Cas A is discussed in Section 3. The current radioactivity 
of ''^Ti in SN 1987A is briefly argued in Section 4. 
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Table 1 

Change of the decay rate and the radioactivity for highly ionized *^^Ti. In the above, 
Ne is the number of bound electrons per atom, and X^^^^ the effective decay rate 
when a ^^Ti isotope has Ng bound electrons. The radioactivity including ionization 

effect, {A + AA)/A, relative to that based on the laboratory decay rate, is calculated 
from equation (7) for both Cas A and SN 1987A. 

2 Activity Change by Ionization: A Linear Analysis 



As mentioned previously, the decay rate of ^^Ti depends on its electric envi- 
ronment. The decay rate. A, is proportional to the inverse of the halflife, 

X='^. (1) 

tl/2 



Given an ionization state, we can compute the electron-capture rate relative 
to the laboratory value as precisely as we like. However, the following approx- 
imation is good enough within the accuracy in question: 

X^Xk + Xli, (2) 



where A is the total decay rate, that is observed in laboratory. The quantities 
Xk and Xli are the partial decay rates capturing K shell (lsi/2) and Lj shell 
(2S1/2) electrons, respectively. The ratio of these partial decay rates for highly 
ionized case is given as 



in a nonrelativistic approximation with a point-charge field. Adopting equa- 
tions (2) and (3), we calculate the effective decay rate for highly ionized ^^Ti. 
This is shown in Table 1 relative to the laboratory value. 

In Table 1, retardation of the decay is manifest in particular when ^"^Ti is in 
H-like and He-like ionization stages. Let us now consider the electron binding 
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energies to see if such a high-ionization is possible. We can simply estimate 
the K- and L-electron binding energies of highly ionized atoms with 

Ee = X 511 [keV] (4) 



again under the assumption that the nucleus is a point-charge and the electrons 
are treated non-relativistically. In the above, a is the fine-structure constant, 
Z the nuclear charge, n the principal quantum number of an electron shell. 
With equation (4) , the binding energies of K shell and Li shell electrons of ^^Ti 
for highly ionized case are calculated to be 6.6 keV and 1.6 keV, respectively. 
It is naturally expected that bound electrons with the range of these binding 
energies can be unbound by shock heating in SNRs seen in X-rays: Even if 
the temperature of a SNR is below E^, the tail of Maxwellian distribution of 
free electron velocity plays a role in ionizing the elements (see Mochizuki et 
al. 1999). 

It should now be clear that the change of ^^Ti radioactivity by high-ionization 
in young SNRs is significant to be investigated. The observable, radioactivity 
A, is generally expressed as 

"^^^TVoAe-^*. (5) 



dt 



Here N is the number of a radioisotope that is synthesized in a supernova 
explosion, Nq is its initial value, and t is the age of a SNR. With an observed 
line fiux E^ the radioactivity is also written as 

A^^^, (6) 
-'7 Ji 



where d is the distance to a SNR, and the absolute intensity of the fiux 
per decay of the parent nucleus. The quantity is the escape fraction of the 
7-photons, that is definitely equal to 1 for Cas A and close to 1 for the later 
phase of SN 1987A. Combining equations (5) and (6), one can easily see that 
the initial mass of ^^Ti is derived from the (observed) values of F^, ti/2, d, and 
the age of the remnant. 

Finally, a linear analysis of the radioactivity (equation [5] ) shows 

^^(l-A.^, (7) 

where AA is the change of the decay rate and AA is that of the activity. It is 
worth noting that AA is always negative, since the ionization always reduces 
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its decay rate. Hence the sign of AA is determined by that of the term in the 
parenthesis in the right side of equation (7). This means the following: If a SNR 
is older than the ^^Ti lifetime (not halflife), ~ 89 yrs, the activity is enhanced 
by ionization, and if younger, the activity is reduced. Very intriguingly, we 
have found that the ionization phenomenon itself affects oppositely in Cas A 
(t = 320 yrs) and SN 1987A (t = 16 yrs). The values of these activity changes 
are also shown in Table 1. We sec in Table 1 that the radioactivity in Cas A 
is enhanced by a factor 2.4 at present and that in SN 1987A is decreased by 
~45% at present if all the ^^Ti atoms are in H-like ionization stage. 

Actual conclusion of the effect of ionization on the radioactivity requires the 
knowledge of the temperature and the density evolution of a SNR. We are 
going to see that the value obtained by the linear analysis is consistent with 
numerical calculations for Cas A in the following session. 



3 ^^Ti Radioactivity in Cas A 

The present radioactivity is affected from the past; the point is whether in 
a SNR ^^Ti could be highly ionized and thus more stable for a considerable 
period of time during the evolution. The basic model adopted here is the 
same as that described in Mochizuki et al. (1999), in which thermal electron 
collisions caused by the reverse shock ionize ^"^Ti as well as ^^Fe. 

We have calculated the ^^Ti radioactivity including the retardation of the 
decay as a function of time (the age of a SNR) and the position in a SNR. The 
calculations are done by solving a hydrodynamical evolution model (McKee 
and Truelovc 1995) with newly introduced clumpy structure, combined with 
microscopic (nuclear/atomic) physics where Dirac-Hartree-Slater method with 
finite nucleus is used to include electron correlations precisely. Our model 
employed here is under update to include so-called di-electric recombination 
process; a free, secondary electron can be bound to a nucleus easier when the 
first electron is already trapped in the orbit. Detailed numerical calculations 
including this process will be presented elsewhere, but inclusion of this process 
will not change the result given below essentially. 

The present result is obtained for the set of parameter values that arc consis- 
tent with recent X-ray observations of Cas A (Willingale et al. 2002): mass of 
the ejecta is taken to be 2Mq and ambient hydrogen density to be 15 cm"^. 
For comparison with a theoretical study of Rauscher et al. (2002), the explo- 
sion energy is set to be 2 x 10^^ ergs for the consistency with their derivation 
of the largest theoretical yield (5 x lO^^Mo) of "^^Ti. The density enhancement 
factor of the clumps (Mochizuki et al. 1999) that contain ^"^Ti is given to be 10, 
which is also compatible with the abundance nonuniformity of the elements 
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(Fe-K and Ni) reported in Willingale et al. (2001). It has been found that rela- 
tively outer region in a SNR is remarkably affected by the ionization and that 
this region almost coincides with the region in Cas A where Fe-K X-ray flux 
is observed from highly ionized Fe. Since ^^Ti is most likely accompanied by 
this Fe, we regard the observed region of Fe-K X-rays (WiUingale et al. 2001) 
as the region where ^'^Ti exists in the remnant. Averaging the radioactivity 
distribution in the remnant over this region, we obtain the averaged enhance- 
ment factor of the radioactivity relative to that without the ionization effect 
as high as ~ 2 at the present age of Cas A. 

It has been frequently argued that theoretical predictions of the initial mass 
of '^''Ti are reasonably smaller than the "observed" initial mass in Cas A as 
inferred from the 7-line measurements on the grounds of the laboratory decay 
rate. Here the "observed" initial mass is only apparent when one does not 
count the ionization effect: The real initial mass is obtained by dividing the 
"observed" mass by the averaged enhancement factor. Accordingly, we see that 
the ionization effect reduces the discrepancy between the "observed" value and 
theoretical predictions. 

We found that the obtained averaged enhancement ratio of the radioactivity, 
~ 2, is large enough to remove the discrepancy in the '^^Ti yield derived from 
the observed values and that in Rauscher et al. (2002) if the real flux value 
is located in the smaller part (< 2.5 x 10~^ photons cm~^s~^) of the reported 
fluxes (Schonfelder et al. 2000 and Vink 2001) within their uncertainties. How- 
ever, the factor cannot compensate the discrepancy sufficiently especially for 
the case that higher half (> 3.5 x 10~^ photons cm~^ s~^) of the reported 7-ray 
flux (Schonfelder et al. 2000) would be real. 

We beheve that the discrepancy is explained largely by the ionization effect. 
However, the remainder of the discrepancy after the subtraction of the ion- 
ization effect, if any, may be attributed to the multi-dimensional effect of the 
explosion: The production of "^"^Ti in 2-D and 3-Dimensional calculations may 
become larger than that in spherical explosion models. The future observa- 
tions of the nuclear fine fluxes and the ionization states of ^^Ti will settle up 
this problem. 



4 ^^Ti Radioactivity in SN 1987A 



We performed Monte-Carlo simulations of Compton degradation of the nu- 
clear 7-photons emitted from the decay sequence of ^^Ti to explain the up- 
per/lower bolomctric luminosity observed at 3600 days after the explosion 
(Suntzeff 1997). Note that at the period of the observation the ionization of 
^^Ti is not relevant. Some details of our calculation are found in Kumagai et 
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Table 2 



Prediction of the nuclear fluxes associated with the Ti decay in 1987A for 6000 
days after the explosion (2003). The ^^Ti mass is given in the unit of Mq, and the 
fluxes are in photons cm-^ s-^. 

al. (1993); nuclear decay parameters adopted in the present study have been 
updated. 

We have obtained the initial ^^Ti mass of (0.82 — 2.3) x 10"^ Mq within known 
accuracy of the experimental values: ti/2 = 60 ± 3 yrs (3(7 deviation) and 

the distance to SN 1987A, d = 48.8 ± 3.3 kpc {3a, Gould & Uza 1998). The 
expected nuclear fluxes for 6000 days after the explosion (i.e., 2003) for both 
the upper and the lower '^'^Ti masses are summarized in Table 2; the escape 
fraction (equation [6]) of the 7-ray photons, depending on each photon energy, 
has been taken into account here. Note that the derived ^^Ti mass depends 
on the distance but the expected fluxes are not. 

In the end, we point out that the ionization process of "^^Ti is considered to 
be well underway in SN 1987A, due to shock heating caused by the collision 
of the supernova blast shock with the dense inner ring. Note that H-like and 
He-like ionization stages of O, Ne, Mg, and Si have been already observed, and 
SN 1987A is a very rapidly evolving remnant (see e.g.. Burrows et al. 2000; 
Michael et al. 2002). If *^^Ti reaches the high- ionization, the expected fluxes 
given in Table 2 become smaller as discussed in the hnear analysis. Details are 
found in Motizuki, Kumagai, and Nomoto (2003). 
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